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COMPUTER  STRATEGY  FOR  DETECTING  LINE  FEATURES 
ON  SIMULATED  BINARY  IMAGE  ARRAYS 
IN  SUPPORT  OF  RADAR  FEATURE  EXTRACTION 


INTRODUCTION 

The  purpose  of  this  study  is  to  develop  meiliods  for  line  feature  detection  and  analysis  on  digi¬ 
tal  n^ar  imagery.  Line  features  on  radar  images  usually  appear  as  a  line  of  constant  gray  tone 
or  as  boundaries  between  area  features.  They  are  characterized  by  length,  width,  shape,  orien- 
tatirni,  and  gray  ume.  Line  features  on  radar  imagery  may  represent,  for  example,  highways, 
railroads,  bridges,  canals,  rivers,  shorelines,  mountain  ridges,  drainage  patterns,  boundaries 
between  agricultural  fields,  and  other  boundaries.  For  automated  line  feature  extraction  it  is 
useful  to  convert  the  radar  image  first  into  a  binary  image.  This  can  be  done  with  image  pro¬ 
cessing  techniques  such  as  edge  detection,  edge  enhancement,  thresholding  and  segmentatioa 
The  binary  image  has  then  to  be  searched  for  line  features.  One  method  to  search  for  line 
features  on  imagery  is  to  search  for  line  paths  in  windows  of  the  image.  An  image  can  be 
divided  into  windows  of  equal  size;  for  example  an  image  consisting  of  512  by  512  pixels  can 
be  divided  into  1024  windows  of  the  size  16  by  16  pixels.  The  line  path  of  a  window  starts  at 
one  side  of  the  window  and  ends  at  another  side  of  the  window.  A  path  of  a  window  is  a 
string  of  pixels  which  arc  connected  by  4-conncctiviiy  except  for  the  end  points  of  the  path. 
The  path  through  the  window  can  be  approximated  by  a  line  clement  characterized  by  its  slope 
and  location  coordinates  in  the  window,  llic  slopes  of  all  line  elements  of  the  image  can  be 
organized  into  a  slope  histogram  or  registered  as  another  image  array.  The  window  line  ele¬ 
ments  can  be  analyzed  for  their  properties  of  forming  lines  of  definite  length,  shape,  and  for 
line  connections  and  parallelism.  This  investigation  deals  primarily  with  search  strategics  for 
finding  paths  in  image  windows.  The  binary  window  arrays  for  the  study  were  computer  gen¬ 
erated  rather  than  derived  from  real  images  because  they  can  be  easily  produced  and  modified 
in  a  controlled  manner. 


LINE  FEATURES  IN  RADAR  IMAGES 

Figures  1  through  6  show  radar  images  that  represent  different  landforms  and  cultural  features, 
including  line  features.  The  radar  images  of  figures  1  through  3  are  from  locations  in  Ala¬ 
bama;  the  images  in  Figures  4  through  6  arc  from  locations  in  Alaska. 

Figure  1  shows  a  IcKaiion  that  contains  mountains  and  a  valley  with  cultural  features.  At  the 
left  side  of  the  image  one  can  recognize  a  dark  line  that  ruas  from  the  lop  to  llic  bottom  of  the 
image.  This  line  feature  is  a  road.  In  the  mountain  area  one  can  detect  at  least  two  major  line 
features  running  from  lop  to  bottom  and  representing  a  valley  between  mountain  ranges.  In  the 
valley  portion  of  the  image  one  can  recognize  many  line  features.  Most  of  the  line  features 
that  arc  short  and  form  right  angles  are  boundaries  between  agricultural  fields.  The  pro¬ 
nounced  bright  line  features  in  the  valley  arc  lines  of  imes  or  shrubs  along  roads  or  creeks  or 
they  arc  boundaries  between  agricultural  fields.  The  tree  lines  that  arc  along  roads  or  fields  arc 
usually  straight.  Those  along  creeks  arc  usually  not  straight  The  line  features  in  Figure  2  are 
produced  primarily  by  valleys,  ridges,  and  roa^!  The  line  features  of  this  image  are  not  ran¬ 
domly  distributed,  but  are  structured  according  to  the  paiicms  of  the  landform. 


In  Figure  3  the  most  obvious  line  features  are  the  two  bridges  across  the  river  and  shorelines 
of  the  river.  The  shorelines  appear  on  a  macroscale  as  straight  lines,  but  show  many  irregulari¬ 
ties  on  a  microscale.  The  city  to  the  right  of  the  river  shows  many  line  features,  which  arc 
caused  by  liiKS  of  buildings  and  roads.  These  line  feauires  have  preferred  directions  in  accor¬ 
dance  with  the  city’s  street  pattern. 

Figure  4  shows  line  features  representing  shorelines,  two  runways  of  an  airfield  appearing  as 
black  lines  in  the  lower  right  comer  of  the  image,  and  straight  forest  boundaries.  The  line 
features  in  Figure  S  arc  produced  by  rivers,  mountain  ridges,  and  shadows.  In  Figure  6  the 
shoreline  represents  a  relatively  straight  line  feature,  and  the  beach  ridges  arc  shown  as  a  fam¬ 
ily  of  gently  curved  line  features. 

The  above  discussed  examples  of  line  features  in  radar  imagery  demonstrate  that  the  lengths, 
orientations,  distribution,  and  geometric  relationships  of  line  features  may  be  used  in  the  pro¬ 
cess  of  detecting,  recognizing,  and  classifying  natural  and  cultural  terrain  features. 

The  fact  that  the  orientation  of  line  features  on  radar  imagery  arc  not  randomly  distributed  but 
show  normally  only  a  few  preferred  orientations  led  to  the  idea  to  develop  line  feature  search 
strategies  for  preferred  orientations. 


INVESTIGATION 

SIMULATION  TESTBED.  A  digitized  radar  image  can  be  partitioned  into  rectangular  win¬ 
dows  that  can  be  represented  by  arrays  of  pixels.  In  this  investigation  the  windows  were 
selected  as  square  arrays,  where  the  number  of  rows  and  columns  is  equal.  The  digitized,  pro¬ 
cessed  image  is  assumed  to  be  binair,  so  that  the  array  elements  are  cither  zero  or  one.  In 
order  to  process  and  manipulate  the  array  elements  a  computer  simulation  testbed  was 
developed. 

The  testbed  algorithm  can  generate  square  arrays  up  to  the  size  of  32  by  32  elements.  The 
arrays  can  be  displayed  or  printed.  The  testbed  provides  the  capability  to  partition  the  array 
into  fields.  The  elements  of  a  field  can  be  selected  so  that  they  arc  all  O’s,  or  I’s,  or  O's  and 
1  ’s  in  a  random  distribution.  A  field  may  also  consist  of  a  line  or  a  path.  The  random  distri¬ 
bution  of  O’s  and  I’s  is  generated  by  a  built-in  random  number  generator.  The  value  of  indivi¬ 
dual  elements  can  be  changed  by  entering  zeros  or  ones  from  the  keyboard. 

With  this  capability  one  can  generate  patlis  consisting  of  ones  that  continue  from  one  side  of 
Uic  window  to  another  side.  All  other  elements  of  the  window  that  do  not  belong  to  an  esta¬ 
blished  path  are  randomly  distributed  zeros  and  ones.  Figure  7  shows  an  array  of  32  by  32 
elements  with  one  path  of  ones  going  from  the  left  side  to  the  right  side  of  the  window.  The 
path  cannot  be  readily  recognized  because  it  is  obscured  by  randomly  distributed  zeros  and 
ones.  In  figure  8  all  elements  except  the  elements  of  the  path  arc  set  to  zeros.  The  path 
represented  by  the  ones  is  easily  recognized. 


SEARCH  STRATEGY.  This  invcsligalion  deals  wiih  the  search  for  paths  that  begin  at  the 
left  side  of  the  array  and  end  at  the  right  side  or  the  bottom  of  the  array.  The  concepts  of  two 
search  codes  and  two  backtrack  codes  arc  developed  and  discussed.  The  search  codes  arc 
called  RIGHT  and  DOWN  aitd  arc  represented  by  the  code  symbols  R  and  D.  The  backtrack 
codes  arc  called  LEFT  and  UP  and  arc  represented  by  the  code  symbols  L  and  U.  The  codes 
R,  D,  L,  U  can  be  combined  and  operated  in  the  combinations  R,  RD,  RDL,  RDLU.  Code 
priorities  and  conditions  to  make  decisions  for  search,  backtrack,  and  switch  ate  developed  and 
discussed. 

The  R  code  searches  for  I’s  in  the  first  column  at  the  left  side  of  the  array  by  moving  from 
clement  (0,0)  to  the  elements  (1,0),  (2,0)  , ...  until  an  element  (j,0)  =  I  is  found.  At  this  point, 
the  R  code  triggered  by  the  condition 

(0.0)=  1) 

switches  to  search  for  I’s  along  row  j  such  that  O.O)  =  1,  0.1)  =  )•  —  until  O.L)  reaches  the 
last  column  on  the  right  side  of  the  array  and  the  path  search  is  completed.  If  R  finds  an  ele¬ 
ment  OW  =  0,  R  returns  to  the  element  0+1. 0)  starts  the  next  search.  The  search  strategy 
of  the  R  code  is  discussed  using  an  example  in  the  following  paragraph. 

Figure  9  shows  a  square  array  of  12  by  12  elements.  The  rows  of  this  array  are  counted  j  =  0, 
1,  2,. ...11  and  the  columns  k  =  0,  1,  2, ....  11.  The  array  contains  a  path  that  begins  on  the  left 
side  at  the  element  (3,0)  and  ends  on  the  right  side  at  the  element  (3,11).  The  R  code  starts  at 
clement  (0,0)  =  0  and  moves  to  element  (3,0)  =  1.  The  search  for  the  tow  of  I’s  starts  at  (3,0) 
and  continues  along  row  j  =  3  until  the  element  (3,11)  =  1  is  reached.  At  this  point  the  path 
search  is  completed.  The  R  search  moves  to  (4,0)  and  searches  for  1  ’s  along  the  k  =  0  column. 
At  (9,0)  =  1,  R  switches  to  search  for  I’s  along  the  j  =  9  row  until  the  element  (9,3)  =  0  is 
found.  At  this  point  the  R  search  of  the  array  is  completed. 

The  RD  code  starts  searching  using  the  R  code.  If  the  clement  (j.k)  =  0  is  found  and  the  con¬ 
dition 

{(j.k)  =  0  and  (j,k-l)  =  1  and  (j+l,k-l)  =  1) 

is  met,  the  R  code  switches  to  the  D  code.  The  D  code  searches  for  I’s  along  the  k-1  column 
such  that  the  conditions 

{(j,k-l)  =  1  and  (j.k)  =  0),  {(j+l,k-l)  =  1  and  (i+l,k)  =  0), 

{a+2,k-l)=  1  and  (j-t-2.k)  =  0). ... 

arc  met  until  an  element  0+d.k)  =  1  is  found.  At  this  point,  triggered  by  the  condition  {(j+d.k) 
=  1),  the  D  code  switches  back  to  the  R  code.  R  and  D  codes  switch  from  one  to  the  other  if 
tiK  required  conditions  arc  met.  The  R  code  has  priority  over  the  D  code.'  Both  search  codes 
arc  combined  into  the  symbol  RD.  The  strategy  of  the  RD  code  is  discussed  with  an  example 
in  the  following  paragraph. 

Figure  10  shows  a  12x12  square  array  similar  to  that  of  Figure  9.  The  array  contains  a  path 
that  begins  on  the  left  side  at  the  clement  (3,0)  and  ends  on  the  right  side  at  the  clement, 
(7,11).  The  search  for  the  path  section  (3,0),  (3,1),  (3,2),  (3,3),  (3,4),  (3,5)  is  executed  by  the 


R  code  and  follows  the  row  number  3.  The  search  for  the  path  section  (4^).  (5^),  (6^).  (7^) 
is  executed  by  the  D  code  and  follows  the  column  number  5.  The  search  of  the  path  section 
(7.6),  (7,7),  (7,8),  (7,9),  (7.10),  (7,11)  is  executed  by  the  R  code  and  follows  the  row  number 
7.  At  (7,1 1)  the  search  for  the  path  is  completed.  The  search  code  RD  begins  with  the  R  code 
and  starts  at  element  (04)).  At  (3.0)  the  condition  {(j.O)  =  1 )  is  met  and  R  switclics  to  a 

search  for  I’s  along  the  row  j  =  3.  At  (3.S)  the  condition 

{(3.6)  =  0  and  (3,5)  =  1  and  (4,5)  =  I ) 

is  met  and  R  switches  to  D.  Tlie  D  code  continues  along  the  k  =  5  column  to  element  (7,5). 
At  this  px)int,  the  condition 

{(7.6)=  1) 

is  met  and  the  D  code  switches  back  to  R  search.  The  R  search  continues  until  the  right  side 
of  the  array  at  element  (7,1 1)  is  reached.  At  this  point  the  path  search  is  completed.  The  RD 
search  moves  to  element  (4.0)  and  starts  the  next  RD  search. 

The  RDL  code  is  used  when  horizontal  dead  end  branches  exist.  The  path  in  Figure  1 1  has 
two  dead-end  bratKhes.  The  elements  (3,6),  (3,7),  (3,8),  (3,9)  form  a  horizontal  dead-end 
branch.  The  elements  (8,5),  (9,5),  (10,5)  form  a  vertical  dead-end  branch.  If  the  RD  code  is 
applied  to  the  array  in  Figure  11,  the  search  will  stop  at  element  (3,10).  In  order  to  find  the 

path,  the  search  code  RD  has  to  be  augmented  by  the  backtrack  code  L.  The  RDL  code 

detects  at  (3,10)  =  0  the  condition 

{(3,10)  =  0  and  (3,9)  =  1  and  (4.9)  =  0) 

which  triggers  the  RDL  code  to  switch  to  the  backtrack  code  L.  In  gciKral  the  RDL  switches 
to  L  if  the  condition 

{a.k)  =  0  and  O.k-1)  =  1  and  O+Lk)  =  0} 

is  met.  The  backtrack  code  L  proceeds  from  element  (j,k)  to  the  elements  (j,k-l),  O.k-2),  ...  if 
the  conditions 

{0,k-l)  =  1  and  O+l.k)  =  0}.  {(j.k-2)  =  1  and  O+l.k-2)  =  0), ... 
arc  met.  If  the  condition 

{0+l.k-l)=  1) 

is  found  the  backtrack  switches  to  the  search  code  D  of  RDL.  In  Figure  1 1  the  L  code  back¬ 
tracks  from  (3,9)  to  (3,5)  where  the  the  condition  {(3,5)  =  1 )  is  found  and  the  switch  from  L 
to  D  occurs.  The  RD  code  searches  until  the  element  (7,1 1)  is  reached  and  the  path  search  is 
completed.  The  second  dead-end  branch  is  ignored  because  R  search  has  priority  over  D 
search. 


4 


The  RDLU  code  is  used  when  a  hori/onlal-vertical  dead  -end  branch  exists  as  shown  in  Figure 
12.  The  array  in  Figure  12  has  two  dead-end  brandies.  The  R  code  of  RDLU  siarts  at  (0,0), 
moves  to  and  continues  along  the  row  j  s  3  to  element  (3,8)  =  0.  The  dead  end  branch 
formed  u;,  the  elements  (1.2)  and  (2.2)  is  ignored  by  the  R  code.  At  (3.8)  the  R  code  switches 
to  tb  D  code.  The  D  code  moves  down  along  the  column  k  =  7  to  the  clement  (6,7)  =  0 
where  the  condition 


{(6.7)  =  0  and  (5.7)  =  1  and  (5.6)  =  0} 

is  found.  This  condition  triggers  the  RDLU  code  to  switch  to  the  U  code.  In  general,  the 
RDLU  code  switches  to  the  U  code  when  the  condition 

{(jjt)  =  0  and  (j-U)  =  1  and  0-l.k-l)  =  0)) 

is  met.  The  backtrack  code  U  proceeds  from  clement  O.k)  to  the  elements  (j-Lk).  0'2.k).  —  if 
the  conditions 


(0-l.k)  =  1  and  (j-l-k-D  =  0).  l0-2.k)  =  1  and  0-2.k-l)  =  0). ... 
arc  met.  If  the  condition 


((j-u.k-l)=  1) 


is  found,  the  backtrack  code  U  switches  to  the  backtrack  L  code.  L  and  U  codes  switch  from 
one  to  another  if  the  required  conditions  arc  met.  In  Figure  12  the  U  code  backtracks  from 
(5,7)  to  (3,7)  where  the  condition  ((3.7)  *  1)  is  found  and  the  switch  from  U  to  L  occurs. 
The  backtrack  code  L  proceeds  from  (3,7)  to  (3,5)  where  the  condition  {(4,5)  =  1}  is  found 
and  L  switches  to  D.  The  RD  code  continues  to  track  the  path  to  the  element  (7,1 1)  where  the 
search  is  completed. 

The  algorithm  represented  by  the  code  RDLU  provides  4  codes.  The  R  code  searches  only  in 
the  right  direction  and  is  trapped  and  terminated  in  right  dead  end  branches.  The  RD  code 
searches  for  right-down  paths,  ignores  branches  in  the  up  and  left  direction,  and  is  trapped  and 
terminated  in  dead-end  branches  .starting  in  the  right  or  down  direction.  The  RDL  command 
searches  only  right-down,  ignores  brarKhes  in  the  up  and  left  direction,  eliminates  dead-end 
branches  extending  in  the  right  direction  and  is  trapped  and  terminates  in  right-down  dead-end 
branches.  The  RDLU  command  eliminates  mo.st  dead-end  branches.  This  search  strategy  is 
further  demonstrated  and  discu.s.scd  with  a  few  examples. 

DISCUSSION.  Figure  13  shows  an  array  of  the  si/c  12  by  12  with  four  paths  staning  at  the 
left  side  of  the  array  and  ending  at  the  element  (9,11)  at  the  right  side  of  the  array.  The  array 
is  called  Original  Array.  The  four  codes  of  RDLU  arc  applied  to  this  array  each  leading  to  a 
different  path.  All  codes  R,  RD,  RDL  and  RDLU  start  at  the  clement  (0,0)  and  search  for  the 
first  element  (j,0)  =1.  At  this  point  the  search  for  the  path  begins. 

The  R  starts  at  (2,0)  =  1.  The  R  code  tests  the  elements  (2,0)  =  1,  (2,1)  =  1,  (2,2)  =  1,  (2,3)  = 
1,  (2,4)  =  0.  Each  time  a  1  is  detected  in  the  Original  Array  a  1  is  entered  at  the  same  loca¬ 
tion  of  the  Transformed  Array.  The  Transformed  Array  has  the  same  size  as  the  Original 
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Array  but  all  elements  are  initially  set  to  zero.  1'lic  K  search  detects  (2,4)  =  0  and  terminates  the  search 
in  the  Original  Array.  The  termination  of  the  search  triggers  the  RULU  algorithm  to  change  all  I’s 
entered  into  the  Transformed  Array  back  to  zeros.  TItc  R  search  continues  to  search  for  the  next  O.O)  = 
1  which  is  found  at  (7.0).  The  R  search  continues  along  row  7  and  terminates  at  (7,9)  =  0.  The  next 
(j,0)  s  1  is  found  at  (8,0).  The  R  search  continues  along  row  8  and  terminates  at  (8.7)  ^  0.  The  fourth 
R  search  starts  at  (9,0)  and  ends  at  (9,11)  where  it  is  completed.  The  path  from  (9,0)  to  (9,11)  is 
entered  into  the  Transformed  Array.  The  search  algorithm  calculates  the  slope  of  the  path.  Figure  14 
shows  the  path  in  the  Transformed  Array  after  the  R  search  of  the  Original  Array  is  completed.  The 
extracted  path  is  a  horizontal  line  with  a  slope  of  zero.  The  path  consists  of  12  elements. 

The  RD  code  starts  at  (2,0).  The  search  terminates  at  (S,8)  =  0.  The  next  RD  search  starts  at  (7,0)  s  l 
and  terminates  at  (7,10)  =  0.  The  third  RD  search  starts  at  (8,0)  =1.  At  (8,7)  the  R  code  switches  to 
D  code.  At  (9,7)  the  D  code  switches  to  the  R  code  and  completes  the  search  at  (9,11).  Figure  15 
shows  the  Transformed  Array  .  If  the  path  is  replaced  by  a  straight  line  from  (8,0)  to  (9,11)  the  slope 
of  this  line  is  0.09.  The  path  consists  of  13  elements. 

The  RDL  code  starts  at  (2,0)  and  terminates  at  (5,8)  of  the  right-down-right  dead-end  branch.  The  next 
RDL  search  starts  at  (7,0).  At  (7,10)  =  0  the  RDL  code  switches  to  the  L  code.  The  L  code  backtracks 
and  removes  the  right  dead'cnd  branch  (73).  (7,9).  At  (7,7)  the  L  code  switches  to  the  RD  code.  The 
RD  code  continues  to  (9.11)  where  the  path  search  is  completed.  Figure  16  shows  the  path  in  the 
Transformed  Array.  The  slope  of  the  line  element  apivoximating  the  path  is  0.16.  The  path  consists  of 
14  elements.  The  backtrack  includes  2  elements. 

The  RDLU  code  starts  at  (2,0).  The  first  right-down-right  dead-end  branch  is  backtracked  by  the  LU 
code  and  removed.  The  second  right  dead-end  branch  is  backtracked  by  the  L  code  and  removed. 
The  RD  code  contineous  to  (9,11)  where  the  path  search  is  terminated.  The  Figure  17  shows  the  path 
in  the  Transformed  Array.  The  slope  of  the  linearized  path  is  0.64.  The  path  consists  of  19  elements. 
The  backtrack  includes  6  elements. 

Figure  18  shows  an  example  of  a  path  having  a  dead  end  with  multiple  branches.  The  RDLU  code  is 
applied  to  this  array  and  the  path  is  shown  in  Figure  19.  The  slope  of  the  linearized  path  is  0.64.  The 
consists  of  19  elements.  The  two  backtrack  branches  include  9  elements. 

The  RDLU  code  represents  only  one  search  strategy.  Another  code  that  also  starts  at  the  left  side  of 
the  array  is  the  DRUL  code.  Ip  this  code  D  has  priority  over  R  and  U  priority  over  L.  Figure  20 
shows  an  array  to  which  the  RD  search  of  the  RDLU  code  and  the  DR  search  of  the  DRUL  code  are 
applied.  Figure  21  shows  the  path  found  with  the  RD  search  and  Figure  22  the  path  found  with  the  DR 
search. 

Figure  2.3  shows  a  priKCSscd,  binary  radar  image  of  the  airport  at  Elizabeth  City,  North  Carolina.  The 
binary  image  con.sists  of  480  x  320  pixels  and  i.s  divided  into  150  windows  each  of  the  size  32  x  32 
pixels.  The  majority  of  windows  have  no  paths  that  extend  from  the  left  side  to  the  right  side  or  bot¬ 
tom  of  the  windows.  The  application  of  the  RDLU  cotic  to  die  windows  produced  23  line  elements 
with  a  siqie  of  0.19  and  12  line  elements  with  a  slope  of  1.15. 

RULES  FOR  SEARCH  CODES.  The  codes  R,  RD,  RDL,  RDLU  start  at  the  element  (0,0)  and  the  ini¬ 
tial  column  (j,0)  of  the  array.  The  search  begins  at  the  elements  O.O)  =  1.  If  the  search  of  any  code  is 
terminated  or  completed  the  code  moves  to  the  next  (j,0)  for  the  next  search. 

The  R  code  searches  for  I’s  along  in  a  row  in  the  right  direction.  If  an  element  (jjc)  =  0  is  detected  the 
search  is  terminated.  The  search  is  completed  when  the  left  column  of  the  array  is  reached.  The  R 
code  ignores  side  branches. 


The  RD  code  combines  the  codes  R  and  D.  R  has  priority  over  D.  The  RD  code  searches  for  I’s 
along  a  row  in  the  right  dirwtioii  and  along  columns  in  the  down  direction.  Specific  conditions  trigger 
the  RD  code  to  switch  from  R  to  D.  from  D  to  R.  to  terminate  the  search,  or  to  complete  the  path 
search.  These  conditions  are  discussed  in  the  section  Search  Suategy.  The  RD  code  ignores  dead  end 
and  side  branches  in  the  left  and  up  direction. 

The  DRL  code  combines  the  codes  R,  D.  and  L.  The  R  and  D  codes  are  the  search  codes  of  RDL,  the 
L  code  is  the  backtrack  code  of  RDL.  The  L  code  backtracks  along  a  row  section  of  a  dead  end  branch 
in  the  left  direction.  Specific  conditions  trigger  the  RDL  code  to  switch  from  R  to  D,  from  D  to  R, 
from  R  to  L,  from  L  to  D,  to  terminate  the  search,  or  to  complete  the  path  search.  These  conditions  are 
discussed  in  the  section  Search  Strategy.  The  RDL  code  ignores  dead  end  branches  in  the  up  direction 
and  side  branches  in  the  left  and  up  direction. 

The  RDLU  code  combines  the  codes  R.  D.  L  and  U.  The  R  and  D  codes  are  the  search  codes  of  RDL, 
the  L  and  U  codes  are  the  backtrack  codes  of  RDLU.  L  has  priority  over  U.  The  U  code  backtracks 
along  a  column  section  of  a  dead  end  branch  in  the  up  direction.  Specific  conditions  trigger  the  RDLU 
code  to  switch  from  R  to  D,  from  D  to  R,  from  R  to  L,  from  L  to  D,  from  L  to  U,  from  U  to  L,  to  ter¬ 
minate  the  search,  or  to  complete  the  path  search.  These  conditions  are  discussed  in  the  secticMi  Search 
Strategy.  The  RDLU  code  ignores  side  branches  in  the  left  and  up  direction. 

The  code  RDLU  can  easily  be  changed  to  other  search  codes  by  permutations  of  the  codes  R,  D,  L,  U 
and  defining  the  initial  column  or  initial  row.  The  search  from  the  top  to  the  bottom  or  to  the  right  side 
of  the  array  can  be  accomplished  by  the  code  DRUL  and  the  initial  row  (OJc).  The  search  from  the  left 
side  to  the  right  side  or  to  the  top  of  the  array  can  be  accomplished  by  the  code  LURD  and  the  initial 
column  (j,n-l),  where  n  x  n  is  the  size  of  the  array.  Code  priorities  and  conditions  for  search,  back¬ 
track  and  switch  are  different  for  each  code. 

SUMMARY,  Cultural  and  natural  terrain  features  frequently  exhibit  line  features  with  preferred  orienta¬ 
tions.  Algorithms  for  directional  search  of  line  features  are  developed  and  discussed.  The  search  sua- 
tegy  includes  two  forward  search  codes  (R  and  D)  and  two  backtrack  codes  (L  and  U).  The  codes  can 
be  combined  to  RD,  RDL  and  RDLU  providing  capabilities  for  determining  parameters  such  as  line  ele¬ 
ment  slope,  path  lengths,  path  end  points  and  number  of  backuack  elements.  These  parameters  arc 
important  to  feature  analysis  and  automated  feathre  extraction. 

CONCLUSIONS 


The  following  conclusions  are  reached: 

Using  directional  search  strategies  a  binary  image  can  be  transformed  into  an  array  of  windows  that 
contain  only  paths  from  one  side  to  another  side  of  the  windows. 

The  path  of  a  window  can  be  approximated  by  a  line  element  that  is  characterized  by  its  slope  and  the 
end  points  on  the  window  sides. 

The  algorithms  for  the  directional  search  suategy  calculate  the  slope  of  the  line  elements,  starting  and 
ending  points  of  the  path,  lengths  of  the  paths  and  dead  end  branches  in  terms  of  number  of  window 
elements. 

The  Transformed  Array  provides  quantitative  information  such  as  length  and  orientation  of  lines,  paral¬ 
lelism,  angles  of  intersecting  lines  and  branching. 
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Figure  1.  Mountains  and  Valley  with  Cultural  Features. 


Figure  3.  City  and  River  with  Bridges. 
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Figure  7.  Window  Array  with  One  Path  Going  from  the  Left  Side  to  the  Right  Side  of 
the  Array.  The  path  u  obscured  by  randomly  distributed  O’s  and  I’s  which  may 
be  representative  for  other  features  or  processing  noise. 
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Figure  8.  The  Same  Window  Array  as  in  Figure  7  with  the  Unobscured  Path.  All  path 
elements  are  ones.  All  other  elements  are  zeros. 
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Figure  9.  Array  with  Right  Path. 
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Figure  iO.  Array  with  Right-Down  Path. 
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Figure  11.  Right-Down  Pith  with  Two  Straight  Dead-End  Branches. 
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Figure  12.  Right-Down  Path  with  Right-Down  Dead-End  Branch. 
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Figure  13.  Original  Array  —  The  Array  Has  Four  Paths  Starting  at  the  Left  Column  and 
Ending  at  the  Right  Column. 
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Figure  14.  Trensformed  Array  —  R  Commend  Applied  to  Original  Array. 
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Figure  IS.  Transformed  Array  —  RD  Command  Applied  to  Original  Array. 
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Fifure  16.  Transformed  Array  —  RDL  Command  Applied  to  Original  Array. 
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Figure  17.  Transformed  Array  —  RDLU  Command  Applied  to  Original  Array. 
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Figure  18.  Original  Array  — 

-  Dead  End  with  Multiple  Branches 
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Figvr*  19.  Transformed  Array  —  RDLU  Command  Applied  to  Original  Array. 
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Figure  20.  Original  Array  Array  with  Paths  Starting  at  the  Left  Side  and  Ending  at  the 
Right  and  Bottom  Side  of  the  Array. 
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1  Figura  21.  TnuMformed  Amy 

—  RD  Commead  Applied  to  Amy  ia  Figure  20. 
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Figure  22.  Trmsfomied  Amy 

—  DR  CoBiauuid  Applied  to  Amy  ia  Figure  20. 
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Figure  23.  Binary  Radar  Image  v>f  Airport  Run  Ways. 
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